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ABSTRACT 

We report on the discovery of HAT-P-12b, a transiting extrasolar planet orbiting the moderately 
bright V w 12.8 K4 dwarf GSC 03033-00706, with a period P = 3.2130598 ± 0.0000021 d, transit 
epoch Tc = 2454419.19556 ±0.00020 (BJD) and transit duration 0.0974 ± 0.0006 d. The host star has 
a mass of 0.73 ± 0.02 M©, radius of 0.70lg'Q^ i?0, effective temperature 4650 ± 60K and mctallicity 
[Fe/H] — —0.29 ±0.05. We find a slight correlation between the observed spectral line bisector spans 
and the radial velocity, so we consider, and rule out, various blend configurations including a blend 
with a background eclipsing binary, and hierarchical triple systems where the eclipsing body is a star 
or a planet. We conclude that a model consisting of a single star with a transiting planet best fits 
the observations, and show that a likely explanation for the apparent correlation is contamination 
from scattered moonlight. Based on this model, the planetary companion has a mass of 0.211 ± 
0.012 Mj, and radius of 0.959j:g:^2i -^J yielding a mean density of 0.295 ± 0.025gcm-3. Comparing 
these observations with recent theoretical models we find that HAT-P-12b is consistent with a ~ 
1 — 4.5 Gyr, mildly irradiated, H/He dominated planet with a core mass Mc < lOAfg. HAT-P-12b 
is thus the least massive H/He dominated gas giant planet found to date. This record was previously 
held by Saturn. 

Subject headings: planetary systems — stars: individual (HAT-P-12, GSC 03033-00706) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

Transiting extrasolar planets (TEPs) provide unique 
opportunities to study the physical properties of plane- 
tary mass objects outside of the Solar System. By com- 
bining time-series photometric observations taken dur- 
ing transit with radial velocity (RV) measurements of 
i the star, it is possible to precisely measure the mass 
and radius of the planet, if the stellar mass and ra- 
dius can be determined by other means. The bulk den- 
\ sity of the planet may then be compared with the pre- 
. dictions of th eoretical planetary structu re mod els (e.g. 
■ i Baraffe et all [2 008: Fo rtnev et all l2007t [Burrows et al.l 
\ l2007t iSeager et al, ,20071 ) to infer the structure of the 
planet. Discoveries of planets that fall outside the pre- 
dicted mass-radius range (e.g. inflated hot Jupiters such 
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as TrES-4b: [Mandushev et al.ll2007D lead in turn to re- 
finements of these models. TEPs also provide unique 
opportunities to study planetary at mosphere s, including 
their composition (e.g. iCharbonne au et al.l [20 021. and 
their thermal profiles (e.g. iKnutson et al.l 12008) . It is 
also possible to measure the projection of the angle be- 
tween the or bital ax is and the stellar spin axis for these 
planets (e.g. iWinn et al. 2005), whic h may be used to 
test theories of planetary migration (jFabrvckv fc Wim] 
[2009h . 

To date more than 40 TEP discoveries have been pub- 
lished, with the majority coming from dedicated pho- 
tometric surveys^^. These planets span a range cover- 
ing more than two orders of magn itude in mass from 
a Super-Earth TEP (Corot-7b; Teeer et all [200ai and 
Supe r -Neptunes (GJ 436b and HAT-P-llb: ICillon et all 
120071 : iBakos et al.l l2009aH to brown dwarf size objects 
(Coro t-3b, X0-3b; iDeleuil et al'][2008l : iJohns-KrulleFan 
l2008f ). Focusing on the low-mass end, we note that the 
three least massive TEPs with well determi ned masses 
(GJ 436b, HAT-P-llb and HD 149026b, iSato et all 
[2005h . are also the three TEPs with the smallest 
radii (excluding Corot-3b) and highest inferred core 
mass fractions. Above this, we begin to see plan- 
ets with a wide range of radii. The planets WASP- 
11/HAT-P-lOb (IBakos et all [20081 iWest et al. 200811,^ 
WASP -6b (.GiUon et al.ll2009D. HAT-P - lb ([Bakos et all 
I2007b|). OGLE-TR-1 1 lb (IPont et all I2004h. WASP- 
15b dWest erall[2?)09h . and XO-2b (|Burke et al.l[200l 
all have radii larger than or c omparable to that 
of Jupiter, whereas HAT-P-3b ([Torres et al.l 120071 : 

http:/ /www. exoplanet.edu/catalog-transit.php 
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iTorres. Winn fc Holmanll2008f) has a radius that is only 
sUghtly larger than that of Saturn. Given the small num- 
ber of TEPs known with M < 0.5Mj, we cannot yet 
say what is the empirical minimum mass of core-less, 
or envelope dominated, gas giant planets. To explore 
the possible transition from envelope dominated to core 
dominated planets, it is necessary to find more low-mass 
TEPs. 

The Hu ngarian-made Aut omated Telescope Network 
(HATNct iBakos et all I200I survey has been a major 
contributor to the discovery of TEPs. Operational since 
2003, it has covered approximately 10% of the Northern 
sky, searching for TEPs around bright stars (8 < / < 
12.5 mag). HATNet operates six wide field instruments: 
four at the Fred Lawrence Whipple Observatory (FLWO) 
in Arizona, and two on the roof of the Submillimeter Ar- 
ray hangar (SMA) of SAO in Hawaii. Since 2006, HAT- 
Net has announced and published 11 TEPs. In this work 
we report on the 12th such discovery. This planet is only 
the fourth sub-Saturn mass TEP announced to date, but 
unlike the other planets, it is of low density, and appears 
to be H/Hc dominated. 

The structure of the paper is as follows. In § [2] we 
summarize the observations, including the photometric 
detection, and follow-up observations. In § [3]we describe 
the analysis of the data, such as the stellar parameter 
determination (§ 13. ip . blend modeling I3.2|) . and global 
modehng of the data (§ 13. 3p . We discuss our findings in 
§111 

2. OBSERVATIONS 
2.1. Photometric detection 
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Fig. 1. — The unbinned light curve of HAT-P-12 including all 
2927 instrumental /-band 5.5-min cadence measurements obtained 
with the HAT-5 (Arizona) telescope of HATNet (see text for de- 
tails), and folded with the period of P = 3.2130598 days (which is 
the result of the fit described in §[3ll. The solid line shows a boxcar 
transit model fit to the light curve. 

The transits of HAT-P-12b were detected with the 
HAT-5 telescope in Arizona. The region around 
GSC 03033-00706, a field internally labeled as 145, was 
observed on a nightly basis between January and July of 
2006, whenever weather conditions permitted. We gath- 
ered 4205 exposures of 5 minutes at a 5.5-minute cadence, 
of which 2927 images were used in the final light curve of 
HAT-P-12. Each image contained approximately 10,000 
stars down to / 14.0. For the brightest stars in the 
field we achieved a per-image photometric precision of 
3 mmag. 

The calibration of the HATNet frames was done uti- 
lizing standard procedures. The calibrated frames were 
then subjected to star detection and astrometry, as 



described in iPal fc BakosI ()2006D . Aperture photome- 
try was performed on each image at the stellar cen- 
troids derived from the 2MASS catalog (jSkrutskie et al.l 
12006( 1 and the individual astrometrical solutions. The 
resulting light curves were decorrelated against trends 
using the External Parameter Decorrelatign tech nique 
in "constant" mode (EPD, see IBakos et aLll2009aD and 
the T rend Filtering Algorithm (TEA, see iKovacs et al.l 
120051 ). The light curves were searched for periodic box- 
like signals using the Box Least Squares method (BLS, 
see Kovacs et al. 2002). We detected a significant sig- 
nal in the light curve of GSC 03033-00706 (also known 
as 2MASS 13573347-1-4329367; a = 13'^57"^33J8^^_5 = 
-h43'i29'"36.7"; J2000; V = 12.84 ± 0.09, IDroeee et al.1 
I2006f l. with a depth of - 20 mmag, and a period of 
P = 3.2131 days. The dip had a relative duration (first 
to last contact) of g « 0.0303 ± 0.0002, corresponding to 
a total duration of Pg w 2.337 ±0.015 hours (see Fig. [J). 

2.2. Reconnaissance Spectroscopy 

All HATNet candidates are subjected to thorough in- 
vestigation before using more precious time on large 
telescopes. One of the important tools for establishing 
whether the transit-like feature in the light curve of a 
candidate is due to astrophysical phenomena other than 
a planet transitin g a st ar is the CfA Digital Speedome- 
ter (DS; iLathaml 11991 ■ mounted on the FLWO 1.5 m 
telescope. This yields high-resolution (R ~ 35, 000) 
spectra with low signal-to-noise ratio sufficient to de- 
rive radial velocities with moderate precision (roughly 
0.5-1 kms~^), and to determine the effective temperature 
and surface gravity of the host star. With this facility we 
are able to reject many types of false positives, such as 
F dwarfs orbited by M dwarfs, grazing eclipsing binaries, 
triple and quadruple star systems, or giant stars where 
the transit signal cannot be due to a planet. 

We obtained 8 observations of HAT-P-12 with the 
DS. The RV measurements of HAT-P-12 showed an rms 
residual of 0.43 km s~^, consistent with no detectable 
RV variation. The spectra were single-lined, showing 
no detectable evidence for more than one star in the 
system. Atmospheric parameters for the star, includ- 
ing the effective temperature Tcff* — 4500 ± 250 K, sur- 
face gravity log 5* = 4.0 ± 0.2 (log cgs), and projected 
rotational velocity v sin i consistent with zero with an 
asymmetri c error of about 1km s~^, were derived as de- 
scribed by iTorres et al.l (j2002l ) . The effective tempera- 
ture and surface gravity correspond to a mid-K dwarf. 
The mean heliocentric radial velocity of HAT-P-12 is 
-40.51 ± 0.21kms^i. 

2.3. High resolution, high S/N spectroscopy 

Given the significant transit detection by HATNet, and 
the positive DS results that exclude obvious false posi- 
tives, we proceeded with the follow-up of this candidate 
by obtaining high-resolution and high S/N spectra to 
characterize the radial velocity variations and to deter- 
mine the stellar parame ters with higher p recision. Using 
the HIRES instrument (jVogt et al.lll994l ) on the Keck I 
telescope located on Mauna Kea, Hawaii, we obtained 22 
exposures with an iodine cell, plus one iodine-free tem- 
plate. The observations were made on 16 nights during 
a number of observing runs between 2007 March 27 and 
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Fig. 2. — Top: Radial-velocity measurements from Keck for 
HAT-P-12, along with an orbital fit, shown as a function of orbital 
phase, using our best fit period (see §[3]l. Zero-phase is defined by 
the transit midpoint. The center-of-mass velocity has been sub- 
tracted. Note that the error bars show the formal errors given by 
the spectrum reduction pipeline and do not include our estimate of 
the stellar jitter. Second panel: Phased residuals after subtracting 
the orbital fit (also see §[3]l. The rms variation of the residuals is 
about 4.82 ms~^. Th ird p anel: Bisector spans (BS) including the 
template spectrum 13.21 1. The mean value has been subtracted. 
Bottom; Relative S values for the Keck spectra. Note the different 
vertical scales of the panels. 

The width of the spectrometer sht used on HIRES was 
0'.'86, resulting in a resolving power of A/AA « 55,000, 
with a wavelength coverage of ~ 3800 — 8000 A. The io- 
dine gas absorption cell was used to superimpose a dense 
forest of I2 lines on the stellar spect rum and es tablish an 
accurate wavelength fiducial (see .Marcv fc Bu tler 1992). 
Relative RVs in the Sol ar System bar ycentric frame were 
derived as described bv lButler et all ( ,1996,) , incorporat- 
ing full modeling of the spatial and temporal variations 
of the instrumental profile. The final RV data and their 
errors are listed in Tab. [1] The folded data, with our 
best fit (sec § [3]) superimposed, are plotted in Fig. [21 

In Fig. [2] we also plot the relative S index. This 
index is computed following the p rescription given by 
IVaughan. Preston fc Wilsonr(|1978D after matching each 
spectrum to a reference spectrum using a transformation 
that includes a wavelength shift and a flux scaling that 
is a polynomial as a function of wavelength. The trans- 
formation is determined on regions of the spectra that 
are not used in computing the S index. Note that the 
relative S index has not been cali brated to the scale of 
IVaughan. Preston fc WilsonI (|1978D . The relative S index 
does not show any significant variation correlated with 
the orbital phase; such a correlation might have indicated 
that the RV variations are due to stellar activity. 

2.4. Photometric follow-up observations 

To confirm the transit signal and obtain high-precision 
light curves for modeling the system we conducted photo- 
metric follow-up observations with the KeplerCam CCD 
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Fig. 3. — Unbinned instrumental i-, z- and g-band transit light 
curves, acquired with KeplerCam at the FLWO 1.2 m telescope on 
the nights of 2007 March 27, 2007 April 25, 2009 February 5 and 
2009 March 6 from top to bottom. Superimposed are the best-fit 
transit model light curves. In the bottom of the figure we show 
the residuals from the fit. Error bars represent the photon and 
background shot-noise, plus the readout noise. 



on the FLWO 1.2 m telescope. We observed four tran- 
sit events of HAT-P-12b on the nights of 2007 March 
27, 2007 April 25, 2009 February 5 and 2009 March 6 
(Fig. 0). On 2007 March 27, 151 frames were acquired 
with a cadence of 90 seconds (75 seconds of exposure 
time) in the Sloan z-band; the observations were inter- 
rupted at mid-transit due to clouds. On 25 April 2007, 
372 frames were acquired with a cadence of 45 seconds 
(30 seconds of exposure time) in the Sloan z-band. On 5 
February 2009, 218 frames were obtained with a cadence 
of 70 seconds (30 seconds exposure time) in z-band. Fi- 
nally, on 6 March 2009, 213 frames were acquired with 
a cadence of 70 seconds (60 seconds exposure time) in 
the Sloan 5-band. This follow-up g-band light curve was 
obtained to further constrain possible blend scenarios 



The reduction of the images was performed as follows. 
After bias and flat-field calibration, we derived an ini- 
tial second-order astrometric transformation between the 
^ 60 brightest s t ars an d the 2MASS catalog, as described 
in lPal fc BakosI (poM ). yielding a residual of - 0.3 - 0.4 
pixels. The primary reason for precise astrometry is to 
minimize the photometric errors that would originate 
from the centroid errors for the individual stars on each 
frame. Aperture photometry was then performed on the 
resulting fixed positions, using a series of apertures. The 
instrumental magnitude transformation was done in two 
steps: first, all magnitude values were transformed to the 
photometric reference frame (selected to be the sharpest 
image), using the individual Poisson noise error estimates 
as weights. In the second step, the magnitude fit was re- 
peated using the mean individual light curve magnitudes 
as reference and the rms of these light curves as weights. 
In both of the magnitude transformations, we excluded 
from the fit the target star itself and the Scr outliers. We 
performed EPD and TFA against trends simultaneously 
with the light curve modeling (for more details, see §[3] 
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TABLE 1 

Relative radial velocity, bisector and activity index measurements of 

HAT-P-12. 
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* The fitted zero-point that is on an arbitrary scale has not been subtracted from 
the velocities. 

^ The values for a^y are the formal uncertainties from the spectrum reduction 
pipeline and do not include our estimate of the stellar jitter. 

This is a relative S i ndex that is sensitive t o variations in S, it has not been 
calibrated to the scale of lVauehan. Preston &: "Wilson 41978i ). 

This is the iodine-free template exposure for which we do not measure the RV 
but do measure the BS and S index. 



TABLE 2 

Photometry follow-up of HAT-P-12 
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0.00073 
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0.00073 


10.51720 i 


187.75095 


-0.00404 


0.00072 
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187.75196 
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Note. — This table is presented in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

^ The out-of-transit level has been subtracted. These magnitudes 
have been subjected to EPD and TFA procedure, carried out si- 
multaneously with the transit fit. 

^ The raw magnitude values without application of the EPD and 
TFA procedures. 



and iBakos et all l2009aD . From the series of apertures, 
for each night, we chose the one yielding the smallest fit 
rms for the light curve. This aperture conveniently fell in 
the middle of the aperture series. The final light curves 
are shown in the upper plots of Fig. [31 with our best fit 
transit light curve models superimposed (see also § ^ . 

3. ANALYSIS 
3.1. Properties of the parent star 



We derived the initial stellar atmospheric parame- 
ters by using the template spectrum obtained with the 
K eck/HIRES instrument. W e used the SME package 
of iValenti &: Piskunovl (Il996f) a long w ith the atomic-line 
database of IValenti fc Fischeil ()2005n . which yielded the 
following initial values and uncertainties (which we have 
conservatively increased to include our estimates of the 
systematic errors): effective temperature Tog^ = 4650 ± 
60K, stellar surface gravity logt;* — 4.75 ± 0.10 (cgs), 
metallicity [Fe/H] = —0.29 ± 0.05 dex, and projected ro- 
tational velocity wsini = 0.5 ± 0.4 km s^^. 

At this stage we could use the effective temperature 
and the surface gravity as a luminosity indicator, and 
determine the stellar parameters based on these two con- 
straints using a set of isochrones. However, the effect 
of log 5* on the spectral line shapes is typically subtle 
and as a result it is generally a rather poor luminosity 
indicator in practice. For planetary transits, the a/R-i, 
normalized semi-major axis and related Pi, mean stellar 
density typicall y impose a stronger constraint on possible 
stellar models fSozzetti "eFaL|[2007l) . The validity of our 
assumption, namely that the adequate physical model 
describing our data is a planetary transit (as opposed to 
a blend) , is shown later in § 13.21 

Using the values of Tcff*, [Fe/H], and log 5^ from the 
SME analysis, and corresponding qu adratic l imb d ark- 
ening coefficients (oz, bz, etc.) from iClareB ( 2004f ). we 
performed a global modeling of the data (? I3.3p . yield- 
ing a full Monte-Carlo distribution of a/i?*. This was 
complemented by a Monte-Carlo distribution of Toffit 
and [Fe/H], obtained by assuming Gaussian uncertainties 
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based on the 1-a error bars of the initial SME analysis. 
For each combination within the large (~ lO"*^) set of 



Fe/H] values , we searched the stellar 
Baraffe et al.) (p97l ITflflSh models for 



a/R^, Toff*, and 
isochrones of the 

the best fit stellar model parameters (such as M^,, R^,, 
logg*, etc). We interpolated these isochrones to the 
SME-based stellar metaUicity of [Fe/H] = -0.29 ± 0.05. 
The majority of the parameter combinations in the 
Monte Carlo search did not match any isochrone. In 
such cases (~ 60% of all trials) we skipped to the next 
randomly drawn parameter set. At the end we derived 
the mean values and uncertainties of the physical param- 
eters based on their a posteriori distribution. We note 
that the spread of the input stellar parameters (based on 
the Gaussian uncertainties) was large compared to what 
the isochrones cover as a function of age, due to the very 
slow evolution of K dwarfs (see Fig. This is partly 
the reason for the 40% match ratio. We also note that 
the match ratio is very sensitive to changing fundamental 
parameters of the isochrones, such as the mixing length 
or the metallicity. 

We then repeated the SME analysis by fixing log g^, to 
the refined value of logg* — 4.61 ± 0.01 based on the 
isochrone search, and only adjusting Teff*, [Fe/H] and 
wsini. This second iteration yielded Tcff* = 4591 ±60K, 
[Fe/H] = -0.36 ± 0.04 and vsini = 1.73 ± 0.5kms-i. 
Curiously, the new Teg* and [Fe/H] values from this sec- 
ond iter ation provide a some what inferior match with 
the B araffe et al.l ([1997, 1998) isochrones, as compared 
to the initial match. Possible reasons for this include i) 
systematic errors in the SME analysis due to the rela- 
tively low SNR of our Keck spectra, ii) increasing un- 
certainty in the SME analysis due to the late stellar 
type of the host star (note that HAT-P-12 has a tem- 
perature that is below th e 4700 K cut-off for sta rs in- 
cluded in the analysis of iValenti fc Fischer! I2005f) . iii) 
general uncertainty in the isochrones for mid-K dwarfs 
(there is a well-known discrepancy between the observed 
and predicted mass-radius relation for K and M dwarf 
stars in double-lined eclipsing binaries such that the ob- 
served radii are larger than the predicted radii, though 
there is some evidence that this discrep ancy only holds 
for rapidly rotating, active stars , e.g. iTorres fc Ribas 
200l iRibai [20061: E6De"z-Moralesl [20071: IChabrier' et al. l 
2007() . Thus we accepted the initial values of Teff*, [Fe/H] 
and V sin i as the final atmospheric parameters for this 
star, along with the isochrone based stellar parameters, 



yielding M^,=0. 733 ± 0.018 Mq, i?*=0.701 



+0.017 
'0.012 -'^O 



Rp, and 



L,t=0.21_Q Q]^ L0. Along with other stellar parameters, 
these are summarized in Tab. [3] 

The stella r evolutionary isochrones from lBaraffe et aTl 
(|1997lll998f) for metallicity [Fe/H]=-0.29 are plotted in 
Fig. m with the final choice of effective temperature Teg* 
and a/Ri, marked, and encircled by the l-cr and 2-cr con- 
fidence ellipsoids. 

The stellar evolution modeling also yields the abso- 
lute magnitudes and colors in various photometric pass- 
bands. We used the apparent magnit udes from the 
2MASS catalogue (jSkrutskie et al.l 120061 ) to determine 
the distance of the system. The magnitudes reported 
in the 2MA SS catalogue have to b e converted to the 
CIT system (|Elias et al.lll982l 119831 ) . in which the stel- 
lar evolution models specify the colors. The reported 




5000 4500 

Effective temperature 



Fig. 4. — Stellar isochrones from lBaraffe et al.l I I1997I . [1999) for 
metallicity [Fe/H]=-0.29 and ages 0.2, 0.5, 1.0, 2.0, 3.0, 4.0, 5.0 
and 6.0 Gyr. The final choice of Teff* and a/R* are marked and 
encircled by the l-cr and 2-cr confidence ellipsoids. 

magnitudes for this star are J2MASS = 10.794 ± 0.023, 
i?2MASs = 10.236 ± 0.022 and X2MASS = 10.108 ± 0.016; 
which are equivalent to Jcit ~ 10.794 ± 0.024, i?ciT = 
10.229 ± 0.023 and Kca t = 10.132 ± 017; in the CIT 
photometric system (see lCarpenteii[200l[ ). Thus, the con- 
verted 2MASS magnitudes yield a color of (J — K) — 
0.66 ± 0.03 that is within l-cr of the expected, isochrone- 
based (J-i^)Baraffc — 0.60 ±0.09. We thus relied on the 
2MASS K apparent magnitude and the Mk = 4.36±0.05 
absolute magnitude derived from the above-mentioned 
modelling to determine the distance: 142.5^3 3 pc. The 
X-band was chosen because it is the longest wavelength 
band-pass with the smallest expected discrepancies due 
to molecular lines in the spectrum of this K4 dwarf. 

3.2. Excluding blend scenarios 
3.2.1. Spectral line-bisector analysis 

Following iTorres et al.l (|2007f ). we explored the possi- 
bility that the measured radial velocities are not real, 
but are instead caused by distortions in the spectral 
line profiles due to contamination from a nearby unre- 
solved eclipsing binary. A bisector analysis based on the 
Keck sp ectra was done as described in §5 of iBakos et al.l 
(|2007ari . 

Fig. [2| shows the bisector spans (BS) phased with the 
orbital period of the planet. In calculating the BS we use 
the convention 

BS = i;iow - I'high (1) 

where viovj is the velocity of the bisector of the cross- 
correlation profile at a low cross-correlation value, and 
Vhigh is the velocity at a high cross-correlation value. 
While the BS do not show significant variations with 
an amplitude that is comparable to or greater than the 
RV variations, there does appear to be a correlation be- 
tween the BS and RV measurements. Applying a Spear- 
man rank-order correlation test, we find that the two 
variables are correlated with 98% confidence (i.e. there 
is a 2% false alarm probability). Since such a correla- 
tion might indicate a blend scenario, we consider below, 
and rule out, the possibility that the system is a blend 
between a bright foreground K star and a background 
eclipsing binary I3.2.2p or a hierarchical triple system 
fij I3.2.3p . The scenarios that we consider are summarized 
in Table m We then consider the possibility that the cor- 
relation is not astrophysical, but rather results from vari- 



6 



Hartman et al. 



TABLE 3 
Stellar parameters for HAT-P-12 



Parameter 


Value 


Source 


Teff* (K) 


4650 ± 60 


SME^' 


[Fe/H] 


-0.29 ±0.05 


SME 


V sin i (km s ^ ) 


0.5 ±0.4 


SME 


Vmac (kms~^) . 


2.26 ±0.0 


SME 


v^ic (kms-i) . 


0.85 ±0.0 


SME 


7rv (kms-i) . 


-40.51 ±0.21 


DS 


az 


0.3432 


SME+Claret^^ 


6z 


0.2493 


SME+Claret 


a,i 


0.4323 


SME+Claret 


bi 


0.2269 


SME+Claret 


^9 


0.8431 


SME+Glaret 


bg 


-0.0064 


SME+Claret 


(Mq) 


0.733 ±0.018 


Baraffe+a/ii;*+SME= 


i?* (i?©) 


701+0 017 


Baraffe+a/iJ*+SME 


log 9* (cgs) .... 


4.61 ±0.01 


Baraffe+a/iJi+SME 


L* (Le) 




Baraffe+a/iJ*+SME 


V (mag) 


12.84 


TASS 


My (mag) .... 


6.89 ±0.11 


Baraffe+a/iJ^+SME 


K (mag,CIT) 


10.132 ±0.017 


2MASS+Carpenter'^ 


Mk (mag,CIT) 


4.36 ±0.05 


Baraffe+a/iJ*+SME 


Age (Gyr) 


2.5 ± 2.0 


Baraffe+a/iJ^+SME 


Distance (pc) . . 




Baraffe+a/iJ*+SME 



observation from the FLWO 1.2 m. We can rule out a 
background object down to ^ 19 mag within ~ 3" of the 
current position of HAT-P-12. Any background object 
must be > 6 mag fainter than HAT-P-12 and thus could 
not be responsible for the observed transit. 



Fig. 5.— Images of a 2'.5 X 2.'5 field containing HAT-P-12 from 
the POSS-I Red survey (Left), and from our FLWO 1.2 m ^-band 
follow-up observations (Right, see section [2.41 1. The dates of the 
exposures are 1955 April 13 and 2009 February 5 respectively. The 
cross marks the position of HAT-P-12 in 1955 and the triangle 
marks the position in 2009. Between these two dates HAT-P-12 
moved ~ 7"5 to the West. From the POSS-I image we can rule out 
the presence of stars brighter than ij ~ 19 at the current position 
of HAT-P-12. 



^ SME = "Spectroscop y Made Easy" packag e for analysis of 
high-resolution spectra I Valenti fc Piskunovl jT996). These 
parameters depend primarily on SME, with a small depen- 
dence on the iterative analysis incorporating the isochrone 
search and global modeling of the data, as described in the 
text. 

^ SME+Claret = Based on the SME analysis and tables of 
quadratic limb darkening coefficients from Claret (2004J. 
" Baraffe+a/i?*+SME = IBaraff'e et al.. ()1997l . il998l ') 
isochrones, a/Ri, relative semi-major axis, and SME re- 
sults. 

Based on the relations from ICarpenteil 11200 II) . 

ations in the sky contamination of the spectra I3.2.4|) : 
we conclude that this is the most likely explanation for 
the correlation. 

3.2.2. Contamination from a background eclipsing binary 

The high proper motion of HAT-P-12 allows us to rule 
out one possible scenario that could potentially fit the 
available observations, namely a background eclipsing 
binary that is aligned, by chance, with the foreground 
K4 dwarf HAT-P-12 (we refer to this as the H,b(S-s) 
model, where H stands for the foreground star HAT- 
P-12, and the b(S-s) is a background eclipsing binary 
star system; here 'b' refers to the fact that the eclips- 
ing binary is in the background rather than being as- 
sociated with the star H). To reproduce the observed 
~ 2.5% deep transit the background object cannot be 
more than 4 mag fainter than HAT-P-12 (objects fainter 
than this would contribute less than 2.5% of the total 
combined light and so could not cause the transit even 
if they were to be completely eclipsed by an object that 
emits no light). Because HAT-P-12 has a high proper 
motion (137.36 ± 2.23 mas yr~ ^ HZacharias et al.ll2004[ ) it 
is possible to use the Palomar Observatory Sky Survey 
plates from 1955 (POSS-I, red and blue plates) to view 
the sky at the current position of H AT-P-12 (this sam e 
technique was used for HAT-P-11. iBakos et al.l [2009al ) . 
Between 1955 and the follow-up observations in 2009, 
HAT-P-12 has moved ~ 7'.'5. Fig. [5] shows an image 
stamp from the POSS-I plate compared with a recent 



3.2.3. Detailed modeling of a hierarchical triple 

Following IBakos et al.l (|2009ar ) we consider the possi- 
bility that HAT-P-12 is a hierarchical triple system, con- 
sisting of two eclipsing bodies that are diluted by a third 
star. In the following we refer to the bright star, with 
properties determined from the SME analysis, as HAT- 
P-12. We consider three scenarios. In the first scenario 
we assume that the bright star HAT-P-12 is uneclipsed, 
and that the two eclipsing components are stars with 
parameters constrained by common origin to fall on the 
same age/metallicity isochrone as HAT-P-12 (we refer 
to this model as the H,S-s model, where H denotes the 
bright star HAT-P-12 and S-s denotes a physically asso- 
ciated eclipsing binary consisting of a brighter star S and 
a fainter star s). In the second scenario we assume that 
HAT-P-12 is uneclipsed, that one of the eclipsing compo- 
nents is a fainter star and that the other component is a 
planet with negligible mass and luminosity compared to 
the star (the H,s-p model, where s stands for the fainter 
star, and p is the planet). In the third scenario we assume 
that HAT-P-12 is a star that is transited by a planet and 
that there is a fainter star diluting the observed transit 
(the H-p,s model). These models will be compared to the 
fiducial model of a single star orbited by a planet (the 
H-p model). 

For the H,S-s and H,s-p scenarios we fit the follow-up 
i-, and g-h&nd light curves together with the HATNet 
/-band light curve following the procedure described by 
IBakos et al.l(|2009aD . We include the HATNet light curve 
to constrain the possibility of a secondary eclipse; to ex- 
clude points that do not contribute to the fit, we only 
include points that are within 1 transit duration of the 
start of transit ingress or end of transit egress, or within 
1 transit duration of the start of secondary ingress or 
end of secondary egress. We use the TFA HATNet light 
curve and apply EPD on the out-of-transit portion of the 
follow-up light curves. We scale the formal photometric 
errors on each light curve so that x^/d.o.f. = 1 for the 
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TABLE 4 

Blend Configurations 



Abbreviation 



Description 



Excluded by 



H,b(S-s) Eclipsing binary star diluted by unresolved, unrelated star 

H,S-s Hierarchical triple star system, two components are eclipsing 

H,s-p Binary star system, fainter star has a transiting planet 

H-p,s Binary star system, brighter star has a transiting planet 

H-p Single star with a transiting planet 



Proper Motion (§ 13.2. 2| l 
Light curve fit (§ 
Light curve fit (§ 
Light curve fit (§ 
Not excluded 



out-of-transit portion of the light curve. We take the 
magnitudes and radii of the stars from the lBaraffe et al.l 
(|1997l 119981 ) isochrones, transforming the BVRI mag- 
nitudes to the Sloan sy s tem u sing the relations from 
iJordi. Grebel fc AmmonI (|2006f ). To make a fair com- 
parison between the blend models and the H-p model we 
also fit the H-p model to the light curves using the same 
procedure used to fit the blend- models (see § 13.31 for a 
more detailed analysis of the H-p model used for the fi- 
nal parameter determinations). Fig. [6] compares the best 
fit H,S-s, H,s-p and H-p models. 

The best fit H,S-s model, consisting of an eclipsing pair 
with masses Mi = 0.73Mq and M2 = Q.12Mq that is di- 
luted by the star HAT-P-12 (M = 0.73 ± 0.02 Mq), has 
= 1445 with 1333 degrees of freedom. We compare 
this to the best fit H-p model which has = 1364 with 
1334 degrees of freedom. Because the photometric noise 
appears to be temporally correlated (Fig. [6|) , formal es- 
timates for the significance of Ax^ between two models 
will overestimate the confidence with which one model 
can be rejected in favor of another. We therefore con- 
duct Monte Carlo simulations to estimate the expected 
distribution of Ax^ values under the assumption that 
the H-S,s model is correct, and accounting for tempo- 
ral correlations in the noise. To generate light curves 
for the Monte Carlo simulations that have similar time- 
correlated noise as the real light curves we Fourier trans- 
form the residual of each light curve from the best-fit H,S- 
s model, randomize the phases, inverse Fourier transform 
it, and then add in the H,S-s model. This method forces 
the simulated light curves to have the same noise power- 
spectrum (and hence auto-correlation function) as the 
actual light curve residuals. We scale the errors of each 
simulated light curve to have x^/d.o.f. = 1 in the out-of- 
transit portion of the light curve. The Fourier transforms 
are carried out assuming uniform time-sampling (a good 
approximation for the follow-up light curves); to use the 
Fast-Fourier Transform algorithm, we cyclically repeat 
each light curve so that the total number of points in 
the light curve is a power of 2. We then fit the H-p and 
H,S-s models to the simulated sets of light curves and 
record Ax^ = Xh p ~ Xh s-s f^'' each simulation. From 
1000 simulations we find a median value of Ax^ = 2.41 
with a standard deviation of 12.4; we find no instances 
where Ax^ < —81, the minimum value attained is —45. 
We conclude that, based on the light curves, the H,S-s 
model can be rejected in favor of the H-p model at the 
> 6(7 confidence level. 

For the H,s-p scenario, we find that the best-fit model 
consists of a star with mass Mi = 0.12Mq transited by a 
planet with Rp = 1.35i?j and diluted by the star HAT-P- 
12 (M = 0.73±0.02 Mq). This model has x^ = 1390 with 



1333 degrees of freedom. To determine the significance of 
— Xh p ~ Xh s-p = ^-^^ '^'^ repeat the Monte Carlo 
simulations, this time adopting the best-fit H,s-p model 
as the fiducial model. From 1000 simulations we find a 
median value of Ax^ = 14.3 with a standard deviation 
of 13.9. There are 4 simulations with Ax^ < —26, so we 
conclude that the H-p model is preferred over the H,s-p 
model with 99.6% (< 3cr) confidence. 

As described in § 13.31 it is possible to correct for sys- 
tematic errors in the photometry by simultaneously ap- 
plying EPD and TFA to the light curves while fitting a 
physical model to them. By using a more sophisticated 
model of this form we are able to rule out the H,s-p 
model with higher confidence, and also rule out the H- 
p,s model. We perform the global modeling as described 
in § 13.31 incorporating three additional parameters that 
allow for dilution in the g-, i- and z-bands. This model 
effectively encompasses both the H-p,s and H,s-p models 
because the only H,s-p models that provide a reasonable 
fit to the light curve are models where the planet-bearing 
star has a mass that is nearly equal to the mass of the 
diluting star HAT-P-12. We allow the dilution factors to 
vary independently in the fit. We find that models with 
no dilution are strongly preferred, and place la upper 
limits on the light contribution in each filter from an un- 
eclipsed star ("third fight") of l^^z < 3%, ^3,^ < 3.5% and 
l3,g < 5%. Any additional star thus makes a negligible 
contribution to the total light of the system. This test 
thus rules out both the H-p,s and H,s-p models. 

3.2.4. Bisector Variations Induced by Sky Contamination 

As shown in sections 13.2.31 and 13.2.21 blend scenarios 
involving an eclipsing binary star system are inconsis- 
tent with other observations of the system. We therefore 
look for an explanation of the apparent BS-RV correla- 
tion shown in § 13.2.11 that does not invoke a blend. One 
possibility is that it is due to varying contamination from 
the sky spectrum. Because HAT-P-12 is relatively faint, 
the flux from the sky is non-negligible compared to the 
flux from the source. Scattered light from the moon il- 
luminating the sky near HAT-P-12 has a solar- like spec- 
trum which yields a peak near RV = Okms~^ in the 
cross-correlation profile of the observed spectrum. The 
degree to which this second peak contaminates the peak 
from the star varies with the sky brightness and the ra- 
dial velocity difference between the moon and the star. 
Because the observer-centric velocity of HAT-P-12 is al- 
ways less than the velocity of the moon, an increase in 
the sky brightness or a decrease in the velocity difference 
will lead to a positive BS variation for our adopted BS 
sign convention. 

While the sky brightness is not directly measureable 
from the available data, in order to quantify this effect we 
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Fig. 6. — Hierarchical triple blend model fits to the follow-up z-band (a) and g-band (b) light curves. EPD filtering was performed on 
the light curves in out-of-transit mode (i.e. a transit model was not simultaneously fit to the light curves). We compare the hierarchical 
triple models to a model consisting of a single star orbited by a planet. The gray-scale points are the un-binned data while the dark points 
are the binned data. Due to time correlations in the noise, the scatter in the binned data is higher than would be expected for white noise 
(the error bars show the expected errors if the noise were white). To save space we only show the 2- and g-band light curves, although 
the follow-up Sloan j-band and HATNet 7-band light curves were also included in the fit. Panels (c) and (d) show the residuals about the 
best-fit H-p model for the follow-up z- and g-band light curves, respectively. Note that the vertical axis for the residual light curves is in 
mmag. Panels (e) and (f) show the distribution of X^_p — Xh s— s ^^'^ ^H-p ~ '^H s-p respectively from 1000 Monte Carlo simulations 
of light curves with correlated noise. The arrows mark the observed values. While the best fit H,S-s and H,s-p models are difficult to 
distinguish by eye from the best fit H-p model, they are statistically rejected in favor of the H-p model at ~ 6cr and ~ 3(t confidence 
respectively. 

may introduce a sky contamination factor (SCF) given 

by 

SCF = ^— — ^ (2) 



Ay2 



where / is the ratio of the flux in the spectrum due to 
the moon to the flux due to the star, is the observer- 
centric radial velocity difference between the moon and 
the star, F = 15.28 kms~^ is the width of the Lorentzian 
function that best fits the mean cross-correlation profile, 
and the form for the denominator is chosen because the 
cross-correlation profile is well fit by a Lorentzian func- 
tion. We estimate / via the relation 



Tjt 



10 



-0.4(B„^Ss) 



(3) 



where T is the total flux received in the region of the spec- 



trum used to compute the BS, t is the exposure time, Tq 
and io are the values for the spectrum with the highest 
count-rate (these are used to account for changes in the 
flux received from the star due to variations in the see- 
ing or transparency) , Bg is the B-band magnitude of the 
star (wc take Bs — 13.84 assuming (B — ]/) ~ 1.0 for a 
dwarf star with Toff — 4650K), and Bm is the effective 
magnitude of the sky due to the moon at the position 
of the star (in an area of ~ 2.5 square arcseconds). To 
estimate B^ we use th e model for the sky brightnes s due 
to moonlight given bv iKrisciunas &: Schaefej ( )1991f ). ex- 
tending it t o the B band b y taking {B — V^ = 0.91 for the 
moon (e.g. 'SchaefeH ll9"98f) and k = 0.19 for the extinc- 
tion coefficient (this is a typical value at the summit of 
Mauna Kea for the Johnson B-band which is roughly 
the region of the spe ctrum used to calculate the BS, 
IKrisciunas et al.l[l987f l. The values for range from 
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18.15 to 19.86. When the moon is below the horizon we 
take Bm = 99.99. 

Fig. [7] compares the SCF to the BS values and to the 
orbital phase. Note that we normalize the SCF to have a 
mean value of 1.0. There is a positive correlation between 
the SCF and the BS. By chance, spectra taken between 
orbital phases 0.5 and 1 had higher sky contamination 
on average than those taken between orbital phases 
and 0.5. When points with SCF > 1.0 are removed, the 
correlation between the remaining BS and RV values is 
no longer significant (the correlation significance is 37%). 
We conclude that this a plausible explanation for the 
apparent BS-RV correlation. 

As a further test on this hypothesis we simulate sky 
contaminated spectra and measure the BS values using 
the same procedure as for the actual spectra. To simulate 
a spectrum we take 

= z{t, Vt-V^) + Ix z{t2, Vt, - -ft,) (4) 

where t is the iodine-free template spectrum of HAT-P- 
12, t2 is the iodine-free template spectrum of HAT-P-13 
scaled to have the same total fiux though the B-band 
as t (HAT-P-13 has T^s = 5638 K, and is thus a bet- 
ter approximation to a solar spectrum than HAT-P-12; 
iBakos et al.l l2Q09bl ) . Vt and Vt, are the barycentric ve- 
locity corrections for the templates, Vi is the barycentric 
velocity correction for spectrum i, is the average ra- 
dial velocity of HAT-P-13, / is given by eq. [H and z{x,y) 
is a function that redshifts the spectrum x by velocity y. 
Figure [H compares the SCF to the BS for the simulated 
spectra. The simulations show a correlation between the 
SCF and BS that is comparable to that seen in fig.[7Ka). 
This confirms that sky contamination may affect the BS 
values at the level that is observed. 

3.3. Global modeling of the HAT-P-12 system 

Our model for the fo llow-up light cu r ves u sed ana- 
lytic formulae based on iMandel fc Agoll (|2002f) for the 
eclipse of a star by a planet, where the stellar flux is 
described by quadratic limb-darkening. The limb dark- 
ening coefficients were derived from the SME re sults 
(§01), using the tables provided bv iClareti (|2004D for 
Z-, i-, and g-bands. The transit shape was parametrized 
by the normalized planetary radius p = Rp/Ri,, the 
square of the impact parameter b^, and the reciprocal 
of the half duration of the transit C/^*- We chose these 
parameters because of their simple geometric meanings 
and the fact that t hese show negligible correlations (see 
IBakos et al.ll2009aD . Our model for th e HATNet data 
was the simplified "P1P3" version of the IMandel &: Agoll 
([2002) analytic functions, for the reasons described in 
[Bakos ct al. (2009a). Following the formalism presented 
by ,Pal (2009) , the RV curve was parametrized by an 
eccentric Keplerian orbit with semi-amplitude K, and 
Lagrangian orbital elements {k,h) = e x (cos w, sinw). 

We assumed that there is a strict periodicity in the 
individual transit times. In practice, we assigned the 
transit number Ntr = to the first high quality follow- 
up light curve gathered on 2007 March 27. The ad- 
justed parameters in the fit were the first transit cen- 
ter observed by HATNet, Tc__ 132, and the last transit 
center observed by the FLWO 1.2m telescope, Tc. +212, 
covering all of our measurements with the HATNet tele- 
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Fig. 7.— (a) BS vs. SCF (eq. EJ. The two variables appear to 
be positively correlated, (b) SCF vs. orbital phase. Spectra taken 
between phases 0.5 and 1 had by chance higher sky contamination 
on average than the spectra taken between phases and 0.5. As 
a result the BS are correlated with the orbital phase, (c) BS vs. 
orbital phase shown separately for points with high and low SCF. 
When points with high SCF are removed, there is no longer an 
apparent correlation between BS and orbital phase. 

scopes, and the FLWO 1.2 m telescope. We prefer us- 
ing Tc, -132 and Tc,+2i2 as adjusted parameters rather 
than the per iod and epo ch for the r eason s discussed by 
IBakos et al.l ([2007cil and iPal et aD (|2008h . The transit 
center times for the intermediate transits were interpo- 
lated using these two epochs and the Ntr transit number 
of the actual event. The model for the RV data con- 
tains the ephemeris information through the Tc,_ 132 and 
?c,+2i2 variables ^al 2009). Altogether, the 11 parame- 
ters describing the physical model were T'c,_i32, ^£,+212, 
Rp/Ri,, b^, C/^*i K, k = ecosw, h = esinw, and three 
additional ones related to the instrumental configuration. 
These are the instrumental blend factor i?inst of HATNet 
which accounts for possible dilution of the transit in the 



10 



Hartman et al. 




1 2 3 4 5 

Sky Contamination Factor 

Fig. 8. — BS vs. SCF (eq. (Jjl for simulated spectra with sky 
contamination. The correlation is similar to that seen for the real 
observations (fig. [7^) . 

HATNet light curve, the HATNet out-of-transit magni- 
tude, Mo,HATNGt, and the relative RV zero-point ^rei- 

We extended our physical model with an instrumen- 
tal model that describes the systematic variations of the 
data. This was done in a simila r fashion to the analysis 
presented in lBakos etlll (|2009aD . Basically, the HATNet 
photometry has been already EPD- and TFA-corrected 
before the global modeling, so we only considered system- 
atic corrections to the follow-up light curves. We chose 
the "ELTG" method, i.e. EPD was performed in "local" 
mode with EPD coefficients defined for each night, and 
TFA was performed in "global" mode using the same set 
of stars and TFA coefficients for all nig hts. The under- 
lying physical model was based on the IMandel fc Agoll 
((2OO2D analytic formulae, as described earlier. The five 
EPD parameters were the hour angle (characterizing a 
monotonic trend that changes linearly over time), the 
square of the hour angle, and the stellar profile pa- 
rameters (equivalent to FWHM, elongation, position an- 
gle) . The exact functional form of the above parameters 
contained 6 coefficients, including the auxiliary out-of- 
transit magnitude of the individual events. The EPD 
parameters were independent for all 4 nights, implying 
24 additional coefficients in the global fit. For the global 
TFA analysis we chose 18 template stars that had good 
quality measurements for all nights and on all frames, 
implying an additional 18 parameters in the fit. We ap- 
ply EPD to the template star light curves using the same 
set of parameters as used for the HAT-P-12 light curves 
before incoporating them in the analysis. Thus, the total 
number of fitted parameters is 11 (physical model) -I- 24 
(local EPD) + 18 (global TFA) = 53, i.e. much smaller 
than the number of data-points (> 1000). 

The j oint fit was performed as described in lBakos et al.l 
(I2009al) . We minimized hi the parameter space by 
using a hybrid algorithm, com bining the d ownhill sim- 
plex method (AMOEBA, see iPress et al.l [i992) with 
the classical linear least squares algorithm. Uncertain- 
ties on the parameters were derived using the Markov 
Chain Monte-Carlo method (MC MC, see H bTd* '200" 
using "Hyperplane-CLLS" chains (jBakos et al . 2009 
The a priori distributions of the parameters for these 
chains were chosen from a generic Gaussian distribu- 
tion, with eigenvalues and eigenvectors derived from the 
Fisher covariance matrix for the best fit value. The 
Fisher covariance matrix is calculated analytically us- 



ing th e partial derivatives given by Ea (l2008l ) and IpH 
(|2009f ). Since the eccentricity of the system appeared as 
insignificant {k = 0.052 ± 0.025, h = 0.007 ± 0.046), 
we repeated the global fit by fixing these to 0. The 
best fit results for the relevant physical parameters are 
summarized in Tab. [H Tab. [5] also lists the RV "jit- 
ter" , which is a component of assumed astrophysical 
noise intrinsic to the star that we add in quadrature 
to the RV measurement uncertainties in order to have 
X^/dof — 1 from the RV data for the global fit. In 
addition, some auxiliary parameters (not listed in the 
table) were: _i32 = 2453792.64889 ± 0.00044 (BJD), 
Tc,+2i2 = 2454897.94147 ±0.00038 (BJD), 7^ei = 32.6 ± 
1.3ms^^(for the Keck RVs, note that this does not cor- 
respond to the true center of mass radial velocity of 
the system, but is only a relative offset). The plane- 
tary parameters and their uncertainties can be derived 
by the direct combination of the a posteriori distribu- 
tions of the light curve, radial velocity and stellar pa- 
rameters. We found that the mass of the planet is 
Mp = 0.211 ± 0.012 Mj = 67.02 ± 3.71 Af®, the radius is 

i?p = 0.959 tgl^a? i?j = lO.lhtofiR® and its density is 
Pp = 0.295 ± 0.025 gcm^'^. The final planetary parame- 
ters are summarized at the bottom of Table [5] 

4. DISCUSSION 

Comparing HAT- P-12b to the theoretical models of 
iBaraffe et all (pOOSi ). we find that the mass and ra- 
dius of the planet are consistent with the 1.0 Gyr, 
Z = Zq = 0.02 non-irradiated model, or with a 1.0- 
5.0 Gyr, 0.10 < Z < 0.50 irradiated model (Fig. ^. 
With an equilibrium temperature of Toq = 963 ± 16 K, 
HAT-P-12b has an equivalent solar semi-major axis of 
Ooquiv — 0.084 AU, so the irradiation received by the 
planet, while not insignificant, is less than what is used 
to calculate the irradiated model (ooquiv = 0.045 AU). 
The inferred metal fraction is exp ected to be cl o ser to 
Z = 0.10 than Z = 0.50 for the IBaraffe et all ()2008f ) 
models if the correct irradiation were used. 

We hav e also compared HAT -P-12b to the theoretical 
models of iFortnev etall (|2007| ). In Fig. [9] we have inter- 
polated these models to agquiv = 0.084 AU, and find that 
the mass and radius of HAT-P-12b are consistent with a 
IOM0 core, 1 Gyr model, and lie between the core-less 
and lOM© core, 4.5 Gyr models. We conclude, therefore, 
that HAT-P-12b is most likely a H/He dominated planet 
with a core of perhaps < lOM®, and a total metal frac- 
tion of < 15%. This makes HAT-P-12b the least massive 
H/He dominated gas giant planet found to date; the pre- 
vious record holder was Saturn. 

It is interesting to compare the properties of HAT-P- 
12b to those of Saturn and HD 149026b, the two plan- 
ets with known radii that have masses closest to that of 
HAT-P-12b. Measurements of the mass of HD 149026b 
range from 0.36Mj to 0.37Mj, while det erminat ions of its 
radius range from .65J?. t to 0.813i?T (ISato et al. 200A 
Charbonneau et al."2006'; 'Torres. Winn & Holman 200a 
Winn et a l. 2008; Nutzman et al. 2009; Carter et al] 



2009D . The planet appears to have a significant core, 
with estimates ra nging from 45M0 to 11AM(^ (see 
iCarter et al.l 120091 . and references therein) , implying a 
high metal fraction of Z > 0.4. Saturn has a mass of 
0.299Mj (,Standishil995f) . equatorial radius of 0.843i?j 
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Fig. 9. — Mass-radius diagram of TEPs (filled squares; blue for those found by HATNet and red for those found by other surveys) and 
solar system planets (triangles). HAT-P-12b is shown as a large filled square on the left. Overlaid ar e th e Baraffe et al. (2008) irradiated 
planetary isochrones for ages of 0.5, 1.0 and 5.0 Gyr and metal fractions of Z = 0.1 and Z = 0.5, the [Baraffe et al. (2008) non-irradiated, 
1.0 Gyr, solar metallicity isochrone, and the Fortney et al. (2O03) 1-0 Gyr and 4.5 Gyr isochrones for planets with core masses of Mc = 
and Mc = lOAfg interpolated to the solar-equivalent semi-major axis of HAT-P-12b. We also show the isodensity lines for 0.4, 0.7, 1.0, 
1.33, 5.5 and 11.9 gcm~^ (dashed lines). HAT-P-12b appears to be well-modeled as a 1.0 — 4.5 Gyr, mildly irradiated planet with a core 
mass of Mq < lOMgj. HAT-P-12b is the lowest-mass H/He dominated gas giant planet found to date. 



(|Seidelinann et all 120071 ). an estimated core mass of 
9Mq < Mc ^ 22M^, a nd a total heavy clement frac- 
tion of 0.14 < Z < 0.29 (jSaumon fc Guill ot 2004). Al- 
though HAT-P-12b is less massive than both HD 149026b 
and Saturn, it has a larger radius than both planets. 
Note that HAT-P-12b does not have a detectable ec- 
centricity, so its large ra dius may not be due to tidal 
heating (in the models bv lJackson. Greenberg fc Barnes! 
120081 however, close-in planets may have tidally inflated 
radii even with eccentricities < 0.01). The large ra- 
dius in comparison with Saturn may be due in part to 
the enhanced irradiation received by HAT-P-12b, and to 
HAT-P-12b potentially having a smaller core mass than 
Saturn. HD 149026b, on the other hand, receives more 
irradiation than HA T-P-12b (aeau i y = 0.025 AU using 
the parameters from [Carter et al.]|2009D . so the differ- 
ence in radii suggests that HAT-P-12b has a substan- 
tially smaller core mass and metal enhancement than 
HD 149026b. 

It is interesting to note that the inferred core mass 
of the three planets appears to correlate with the host 
star metallicity (HAT-P-12 has [Fe/H] = -0.29 ± 0.05, 
the Sun has [Fe/ H] = 0, and HP 149026 has [Fe/H] = 
0.36 ± 0.05 from lSato etaIll2005D. This correl ation has 
been previously not ed by iGuillot et al.l ()2006f ) and by 
I Burrows et al.l (|2007f l , and is perhaps suggestive evidence 



for the core accre tion model of planet formation (e.g. 
lAlibert et al.l[2005l and references therein). Further dis- 
coveries of TEPs with masses comparable to or less than 
that of Saturn are needed to determine whether or not 
this correlation holds. Note from Fig. [5] that the radii 
of planets in this mass regime are more sensitive to the 
core mass than are the radii of more massive planets for 
which a given core mass is a smaller fraction of the total 
planet mass. 

Finally, one might wonder why other planets like HAT- 
P-12b have not been found to date (see Fig. The 
significant ~ 2.5% transit depth of HAT-P-12b is weU 
within the range that is easily detectable for many tran- 
sit surveys, and the 35.8±1.9kms~^RV semi-amplitude, 
while small, is still easily measured with high-precision 
RV spectrometers (though more observations may be 
needed for a robust confirmation, which may slow the 
rate at which these planets are announced). We con- 
clude that of hot gaseous planets with radii similar to 
Jupiter, only a small fraction have masses similar to Sat- 
urn such as HAT-P-12b; the majority have masses sim- 
ilar to Jupiter. With the discovery of HAT-P-12b, we 
estimate that the fraction is ~ 2%, with considerable 
uncertainty. 

We would like to thank the referee, Peter McCul- 
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TABLE 5 

Orbital and planetary parameters 



Parameter 



Value 



Light curve parameters 

P (days) 3.2130598 ± 0.0000021 

Tc (BJD) 2454419.19556 ±0.00020 

Ti4 (days) 0.0974 ± 0.0006 

Ti2 = T34 (days) 0.0125 ± 0.0005 

a/_R* 11.77l°-^5 

C/iJ* 23.57 ±0.12 

Rp/R^, 0.1406 ±0.0013 

0.044iH'2^ 

b = acosi/R-i, 0.21ll!5]|{7g 

i (deg) 89.0 ± 0.4 

RV parameters 

K (ms-i) 35.8 ± 1.9 

fcRv'' 

hnv"^ 

e 

RV jitter (ms-^) 4.82 

Planetary parameters 

Mp (Mj) 0.211 ±0.012 

Rp (Rj) 0.959l» o^9 

C{Mp,Rp) 0.63 

Pp (gcm-3) 0.295 ±0.025 

a (AU) 0.0384 ± 0.0003 

loggp (cgs) 2.75 ±0.03 

Tcq (K) 963 ± 16 

0.023 ±0.001 

(F) (10i2ergs-icm-2) o 1.91 ±0.12 
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T14: total transit duration, time between first and last con- 
tact; T12 = T34: ingress/egress time, time between first and 
second, or third and fourth contact. 
Fixed to 0. 

Correlation coefficient between the planetary mass Mp and 
radius Rp. 

The Safronov number is given by = |(V'esc/Vorb)^ = 
(a/Rp)(Mp/Mtc) fsee lHansen fc Barmanll2007h . 
" Incoming flux per unit surface area. 
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